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Abstract 

An ab initio calculation of the 4 He(e, e'p) 3 H longitudinal response is presented. The use of the 
integral transform method with a Lorentz kernel has allowed to take into account the full four- 
body final state interaction (FSI). The semirealistic nucleon-nucleon potential MTI-III and the 
Coulomb force are the only ingredients of the calculation. The reliability of the direct knock-out 
hypothesis is discussed both in parallel and in non parallel kinematics. In the former case it is 
found that lower missing momenta and higher momentum transfers are preferable to minimize 
effects beyond the plane wave impulse approximation (PWIA). Also for non parallel kinematics 
the role of antisymmetrization and final state interaction become very important with increasing 
missing momentum, raising doubts about the possibility of extracting momentum distributions 
and spectroscopic factors. The comparison with experimental results in parallel kinematics, where 
the Rosenbluth separation has been possible, is discussed. 

PACS numbers: 21.45. +v, 25.10.+S, 25.30.Fj, 27.10. +h 
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I. INTRODUCTION 



Numerous experimental as well as theoretical investigations of (e, e'p) exclusive reactions, 
in both light and heavy nuclei, have been performed extensively in thepast, with the aim 



to extract information about the structure of these systems 
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Q, 19, Q . In particular one has tried to access ground state 
properties of the target nucleus like spectroscopic factors, shell momentum distributions etc. 
However, it is well known that such quantities can only be obtained under the hypothesis 
that the reaction mechanism is dominated by a direct knock-out of the proton and neglecting 
the interaction in the final state. Such assumptions are usually considered more and more 
plausible as the momentum transferred by the electron to the system increases and allows 
to probe "single nucleon" physics. In many nuclei the experimental one-body knock-out 
spectra indeed show very nicely pronounced peaks, hinting to an independent motion of the 
nucleons in such systems. In these cases shell momentum distributions and spectroscopic 
factors have been extracted trying to estimate FSI effects in various ways. Spectroscopic 
factors which are found smaller than 1 (and they are often found considerably far from 1) 
are interpreted as due to large "correlation effects" induced by the residual interaction in 
the ground state of the system. 

Unfortunately up to recently the two fundamental assumptions mentioned above could 
not be checked, because of the impossibility to solve the man-body problem in a quantum 
mechanical consistent way both for ground and continuum states. The recent progress made 
in few-body physics allows us to start investigating these assumptions. For A = 2 and 3 the 



calculations are fully under control 



22| both for ground and continuum states, and the 



problem has been investigated. However, the features of such systems are often considered 
too different from those of a typical "many-body" nucleus, to be taken as testing grounds 
for validating assumptions on heavier systems. 

In the last decade it has been demonstrated that the procedure to calculate reactions 
with the help of integral transforms originally proposed in 1985 j3] can be successfully 
applied in order to overcome the longstanding stumbling block which prevented ab initio 
calculations of high energy reactions involving four nucleons and more 24, 25, 26, 27, 
This has been possible thanks to the integral transform with the Lorentz kernel (we will 
denote it by LIT) proposed in Ref. |2i| and recently applied also to the two-body break-up 



2 



of the four-body system in Ref . |30| . Thus it is possible now to treat the full dynamics of a 
reaction to continuum in a nucleus whose features (binding energy and density) are certainly 
much closer to those of heavier systems than deuteron and triton or 3 He. 

It is the purpose of this work to perform a model study of the role of the full treatment 
of the interaction in the four-body dynamics in the 4 He(e, e'p) 3 H reaction and to discuss 
the plausibility of the direct knock-out and plane wave assumptions. To this aim we use the 



semirealistic potential MTI III 



3 it ] and concentrate on the longitudinal response function, 



where meson exchange currents effects are negligible. This response is accessible experimen- 
tally if one performs a Rosenbluth separation in parallel kinematics, and has been measured 
in a number of experiments 3, 17]. Though the use of a semirealistic potential does not 



allow to draw precise conclusions we believe that a comparison with data may be instructive 
and we will also comment on that. 

In Sec. II the expression for the (e, e'p) cross section is recalled and the formalism describ- 
ing the integral transform approach with a Lorentz kernel to exclusive reactions is reviewed. 
Results are given in Sec. Ill while conclusions are summarized in Sec. IV. 



II. GENERAL FORMALISM 



A. Cross Section 



The sixfold electrodisintegration cross section of 4 He into the two fragments p and 3 H is 
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given by 

d 6 a p j 



33] 



cm [V L S L (uj,q,p p ,9p) + V T S T (uj,q,p p ,9 p ) + VltSlt(u, q,p p , 9 P ) cos. 



dE'dQ e /dp p 

+V TT S TT (uj,q,p p ,9 p ) cos20] . (1) 

Here E' and Q e > denote energy and solid angle of the electron after the reaction, <jm is 
the Mott cross section and (f) denotes the angle between the electron and ejectile planes. 
Energy and momentum transferred from the electron to the nuclear system are denoted as 
oj and q = gq. The quantity p p = (p p , Q p ) denotes the momentum of the proton detected 
in coincidence with the electron, and 9 P is the angle between the outgoing proton and q. 
The Vp are kinematical coefficients and the nuclear dynamics is contained by the structure 
functions Sp. 
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Integration over p p leads to the fivefold cross section 



d 5 a p ,t 



dE'dtt e ,dtt r 



d 6 



a 



P,t 



PI 



dE'd^l e 'dpp 



dE„ 



dp p 



dE n 



Pp 



dE„ 



[V l F l (lj, q, 9 p ) + V T F T {u, q, 9 P ) 



+V lt F lt (uj, q, Op) cos 4> + VttF T t(uj, q, 9 P ) cos 2c, 



(2) 



where E m = uj — T p — T t represents the missing energy (T p and T t being the proton and 
triton kinetic energies). The new structure functions Fp are simply given by 



Fp(uj,q,9 p ) = J Sf3(uj,q,p p ,9 p )dE m . 



(3) 



Notice that, since S/3(uj,q,p p ,9 p ) include the energy conserving (^-function, the integration 
over E m fixes a unique value of p p for each combination of u, q and 9 P . 

The total contribution of the p, t disintegration channel to the inclusive cross section is 



d 3 a p ,t 
dE'dVt pl 



[V L R p /(co,q)+V T R%\co,q 



with 



r 2irp 2 

R p S(u, q) = -Fp{u, q, 9 p ) sin 9 p d9 p . 

J dE m 



(4) 



(5) 



In what follows we concentrate on the longitudinal response Fi(q, U),9 P ), representing the 
response of the system to the electron-nuclear charge interaction. This can be written as 



F L {q,u,6 p )= K*ftt(^t)|p(q)|*«>r, 

M t ,M p 



(6) 



where the four-body ground state is denoted by ty a , and is the continuum final-state 
of the minus type pertaining to the proton-triton channel [34] with the relative proton- 
triton momentum k = kk.. The sum goes over the projections M t and M p of the fragment 
total angular momenta in the final state. The continuum states \l/~ t are normalized to 
S(k — k')SM t M^M p M' ■ The quantity E p t is the final state intrinsic energy 



k 2 

En t = h Ef , 

p,t ^ -r t, 



(7) 



where \i is the reduced mass of the proton-triton system and E t denotes the 3 H ground state 
energy. 



The initial and final states are connected by the nuclear charge operator p which we take 
in its non relativistic form 



4 1 4- r 3 

3=1 



exp (iq ■ Tj 



(8) 



Here r| denotes the third component of the j-th nucleon isospin, rj represents the position 
of the j-th nucleon with respect to the center of mass of the four-body system and G P E is 
the proton electric form factor. In comparing our results with experimental data we will use 
the proton form factor G P E = G P E /(1 + (q 2 — cu 2 ) / Am 2 ) 1 ^ 2 (containing first order relativistic 
correction) with G P E in the usual dipole parametrization. 

The main difficulty in the calculation of Fl is represented by the continuum wave function 
^p t (E Ptt ) in the transition matrix element 

T Pit (E Ptt ) = (*- t (E p , t ) |p| V a ) . (9) 

With the integral transform method [3| with the Lorentz kernel one is able to 

perform an ab initio calculation of this transition matrix element in a large energy range 
without dealing with the continuum solutions of the four-body Schrodinger equation. How 
this is possible has been described in Ref. |23[ and will be briefly summarized in the next 



subsection. Further details can be found in 130, 35, 36 



36, 2l|- 



B. The LIT Method for Exclusive Reactions 



The LIT approach to exclusive reactions consists in calculating transition matrix element 
of the perturbation O between the initial (^/q) an d final (^fj) states 



T f {E f ) = (9j(E f ) 



O 



^0 



(10) 



without calculating tyj(Ef). 

In general denoting with a and b the two fragments containing n a and = A — n a 
nucleons, respectively and with H the full nuclear Hamiltonian, we have the following formal 



expression for ^ f =a b (Ef =a ^ b ) in terms of the "channel state" <fif =ab (Ef =a ^) 

1 



E a)b -is- H 



\<p- b (E a , 



(11) 
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where A is an antisymmetrization operator. In case that at least one of the fragments is 
chargeless the channel wave function (p~ b (E a i6 ) is the product of the internal wave functions 
of the fragments and of their relative free motion. Correspondingly, V in Eq. (jllj) is the sum 
of all interactions between particles belonging to different fragments. If both fragments are 
charged, like in our case, <p~ b (E a ^) is chosen to account for the average Coulomb interac- 
tion between them, and the plane wave describing their relative motion is replaced by the 
Coulomb function of the minus type. Correspondingly, V in Eq. (fTTj) is the sum of all inter- 
actions between particles belonging to different fragments after subtraction of the average 
Coulomb interaction, already considered via the Coulomb function. We write 4>~ b {E a ^) in 
the partial wave expansion form 

oo I 



( E a,b 



$ a (l,...,n a )$ 6 (n a + 1,...,A) 
(2tt) 3 / 2 



i=0 m=-l 



w e (k;r) 
kr 



(12) 



Here $ a (l, n a ) and $&(n a + 1,...,A) are the internal wave functions of the fragments, 
r = (r, Q r ) = R" m — H b cm represents the distance between them, and the energy of the 
relative motion is k 2 /2fi = E a ^ — E a — Eb, where E a and Ei, are the fragment ground state 
energies. The functions we(k;r) are the regular Coulomb wave functions of order £, and 
5g(k) are the Coulomb phase shifts 34j . The internal wave functions of the fragments are 
assumed to be antisymmetrized and normalized to unity, so that the properly normalized 
continuum wave function in Eq. (jllj) is obtained via application of the antisymmetrization 
operator. For n a = 1 this has the form 

A 

i=2 



A 



(13) 



where Vij are particle permutation operators 

When one inserts Eq. (fTTj) into Eq. (fTUj) the transition matrix element becomes the sum 
of two pieces, a Born term, 



T: 



Born 
a.b 



and a FSI dependent term, 



VA 



A O 



1 



(14) 



E a h + ie — H 



O 



(15) 
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While the Born term is rather simple to deal with, the determination of the FSI dependent 
matrix element is rather complicated. Within the LIT approach this term is treated as 
outlined in the following. 

In Eq. (115)) one inserts the completeness relation of the Hamiltonian eigenstates \fr v (E v ) 
(labelled by channel quantum numbers v and normalized as ^ v \^ v i) = 5(u — v')) 

1 



dv(<j>-AE a>h )\VA\y v {E v )) 



Enb + ie — E u 



(* v (E v )\0\*o) . (16) 



T a,b ( E a,b) = 

Defining F a ^ b (E) as 

F a , b (E) = j£dis(<f ) - b (E a , b ) VA m v {E v )){^ v {E v ) O * ) 5(E - E v ) 



(17) 



one has 

T a T(E a 



Fa,b( E ) 



-dE 



iirF a , b (E aib ) + V f 
Je: 



Fg, b (E) 

E„,_h — E 



dE 



;is) 



E7, Ea,b + ie — E 

where E t h is the lowest excitation energy in the system i.e. the break-up threshold energy. 

The function F ab contains information on all the eigenstates ^> v for the whole eigenvalue 
spectrum of H . It is obtained by its Lorentz integral transform 

Fa,b(E) 



L [F a , b ] (a) 



E - (E-a R y + aj 



dE 



(19) 



where 



*iW = {H- a)- x O |* ) , M°) = {H- a)- l AV\<l>lb{E a ,b)) (20) 
and a = a R + icrj. Equation (jTUjl shows that L [F a>b ] (a) can be calculated without explicit 



knowledge of F ab , provided that one solves the two equations 



(H-a) 
(H-a) 



$2 



O|* ) , 
AV\4>l h {E a , b )) 



(21) 
(22) 



which differ in the source terms only. It is easy to show that $1 and \&2 have finite norms. 
When solving Eqs. (|2T|) and (|22j) it is sufficient to require that the solutions are localized, 
and no other boundary conditions are to be imposed. Therefore "bound state" techniques 
can be applied. 

We use an expansion over a basis set of localized functions consisting of correlated hy- 
perspherical harmonics (CHH) multiplied by hyperradial functions. As discussed in (26 1 for 



the case of the total 4 He photoabsorption cross section, special attention has to be paid to 
the convergence of such expansions. A rather large number of basis states is necessary in 
order to reach convergence, thus leading to large Hamiltonian matrices. Instead of using a 
time consuming inversion method we directly evaluate the scalar products in (|19|) with the 



Lanczos technique as explained in Ref. 



38|. 



After having calculated L[F a ^](u) one obtains the function F a ^(E), and thus T ai b(E a} b), 
via the inversion of the LIT, as described in 3^|. 

In the next section results obtained by means of Eq. ()14|) will be labelled by PWIAS. 
The label PWIA will indicate that in Eq. ()T4*)l the antisymmetrization operator A has been 
neglected. We remind the reader that in this case the structure function Fff turns out to 
be factorized in terms of the proton form factor and a function n(|q — p p |), which is the 
Fourier transform of the overlap between the 4 He and 3 H ground state wave functions. 



III. RESULTS 



As already mentioned, the ground states of 4 He and 3 He as well as the \& in Eqs. (J21)) and 

n order to speed up the convergence, 



We use the MTI-III 



potential 



(|22j) are calculated using the CHH expansion methoc . 
state independent correlations are introduced as in 2l 
and identical CHH expansions for the ground state wave functions of 4 He and of the three- 
nucleon systems as in Q and Q, 

respectively. 

We calculate the transition matrix elements (j!4|) and ()15|) in the form of partial wave 
expansions. When one substitutes the expansion (JT2J) and the expansion 

p( q ) = ^r L V^)p LM (g) (23) 

LM 

of the charge operator (JHJ) into the Born amplitude (|T4jl and into the right-hand sides of 
Eqs. (12 lj) and (|2~2"|) one finds that in our case of central NN forces the transition matrix 
element (JHJ) turns into a sum over over L (L is equal to the I in 1)12)1) of partial transition 
matrix elements multiplied by the factors 

L 

Y LM {n k )Y£ M {n q ) = (4tt)- 1 (2L + l)P L (k • q). (24) 

M=-L 

These factors determine the dependence of the cross section on 9 p . The dynamic equations 
are split with respect to orbital momentum L and they are M-independent. The multipole 
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transitions of the charge operator in Eq. (JHJ) are taken into account up to a maximal value 
of Lsom = 20 and Lpsi = 6 for the Born and FSI terms, respectively. (The relatively low 
value of Lpsi is due to the fact that FSI does not affect the final state higher partial waves 
significantly.) Correspondingly Eqs. (J21j) and (J22|) are solved for the different values of L, 
running from to Lpsi- Since the excitation operator induces both isoscalar and isovector 
transitions, the hyperspherical harmonics (HH) entering the calculation are characterized 
by the quantum numbers L = 0, 5* = and T = 0, 1. In the calculation involving L up to 4 
the maximal value of the grand-angular quantum number K max is taken 7 (odd multipoles) 
or 8 (even multipoles), the only exception being the L = 1 multipole in the T = channel, 
for which K max = 9 has been used. For L = 5 and L = 6 K max is taken equal to 9 and 10, 
respectively. These values of the grand-angular quantum number provide the convergence 
of the various LIT's of Eq. (|T9"j) with an uncertainty in the response function © of less than 
1%. In addition for K max = 9 and 10 a selection of states has been performed with respect to 
the permutational symmetry types of the HH. Among the HH entering, the expansion, those 
belonging to the irreducible representations [f] = [2] and [f] = [-] j3, 41 1 of the four-particle 
permutation group S4 can be neglected in the calculation of the LIT for K values higher 
than 7 (odd multipoles) and 8 (even multipoles). 

We start illustrating the contributions of the proton-triton channel and of the mirror 
channel due to the neutron- 3 He break-up to the total inclusive response function. This 
comparison serves as a test of the correctness of the results. In fact below the threshold 
for the disintegration of 4 He into proton, neutron and deuteron for the isovector case and 
into two deuterons for the isocalar case, the two results should coincide. The neutron- 3 He 
response R n £ h {uo,q) has been calculated along the same lines described above, except that 
in the "channel state" 4>J =a h {Ef =a ^) of Eq. (|TT|) the relative motion is given by a plane 
wave. We choose to compare the sum of R^iuj^q) and K1' h (uj,q) for the multipoles L=0, 
T=l and L=2, T=0 (two of the multipoles which contribute most) with the total inclusive 
response calculated for the same multipoles. In Fig.^this comparison is shown. Considering 
that the calculation of the total longitudinal response proceeds in a very different way, i.e. 



only by inversion of the norm of \&i 



24j, this comparison confirms the correctness of the 



calculation. Besides the degree of accuracy of the results one notices that for these multipoles 
the proton-triton and neutron- 3 He channels dominate much beyond those thresholds. 
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A. Parallel kinematics 



Our study of Fi(uj,q,9 p ) focuses first on the parallel kinematics of Ref. |l7[ where a 
Rosenbluth separation has been performed. In the columns 2-4 of Table |I] we list the values 
of q, u and modulus of missing momentum p m = q — p p of the kinematics we have chosen 
to analyze (labelled by Kin. N. in column 1). The values of the experimental energies and 
momentum transfers are illustrated in Fig. 2 as points in the q — to plane and labelled with 
the corresponding numbers. In the same figure their positions with respect to the ridge 
u = q 2 /(2m p ) are shown (m p is the proton mass). The value of the final state intrinsic 
energy E Pjt , which is the input of the calculation, has been obtained by calculating first the 
relative momentum k from the relation 

k=J^- {q -^ ) (25) 

and then using Eq. (JJJ). 

In column 5 of Table |l] the PWIA results are listed. In this approximation and in an 
independent particle model of 4 He the PWIA result represents the probability that the 
proton in the S-shell of 4 He has momentum p m . Therefore one has constant values for Kin. 
N. 1-3 and 4-8. The integral over all values of p m gives the "spectroscopic factor" for that 
shell, which for this potential turns out to be 0.88 (3] (this value can be compared with 
0.84 obtained using a realistic potential like AV18 and Urbana IX 44 , 4s|). 

In Table H] the effects of antisymmetrization and of FSI are also shown as percentages of 
the PWIA values (the results denoted as FULL include both effects). In general one notices 
small effects of antisymmetrization for almost all cases as one would expect for kinematics 
with p m much smaller than q. Nevertheless for the kinematics at lower energies these effects 
can increase up to about 10%. The role of FSI is much more important, especially at low q. 
One notices that i) for the kinematics close to the u = q 2 / (2m p ) ridge FSI effects decrease 
for increasing q; ii) Kin. N. 4 and 9, which are more distant from the ridge u = q 2 /(2m p ), 
present a rather high contribution of FSI; iii) at higher momenta and in the lower energy 
side of the ridge FSI enhances the PWIA results. This effect goes in the opposite direction 
compared to previous estimates based either on optical potentials and orthogonalization 
procedure 4f|, or on diagrammatic expansions j^. 

The observation iii) is consistent with previous ab initio calculations of the inclusive 
longitudinal response function in 2 H [21], 3 He 0]) and 4 He j^J. In the latter case one 
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finds that the longitudinal responses at constant q values, calculated with and without FSI, 
cross at an u value of approximately q 2 / (2m p ). The fact that the crossing happens just along 
that ridge is probably due to the different effects of the potential in the initial and in the 
final state with respect to the free one-body knock-out model, as explained in the following. 
In the one body knock-out model the PWIA peak energy is uj pea k = q 2 /(2m p ) + A. The 
positive quantity A is the difference between the binding energies of 4 He and 3 H and can 
be considered as a "ground state effect" of the potential. One can argue that the additional 
effect of the potential in the final state would lead to uj peak = q 2 /(2m p ) + A — V, where V 
represents the mean interaction energy between the proton and the triton in the final state 
interaction zone (V will be attractive). Therefore the PWIA and FSI curves should intersect 
at an energy smaller than q 2 / (2m p ) + A. To a good accuracy this value turns out to be just 
q 2 /(2m p ). Of course such a comparison between inclusive and exclusive results is justified 
only in case of sufficiently low p m as it is the case for the kinematics listed in Table |1J 

Similar PWIAS and FSI effects are also found for Kin. N. 6 and 9 in the two-body 
break-up results of 3 He 

It is a common belief that the kinematical regions at lower energy and higher momentum 
transfers are the privileged ones to investigate the ground state short range correlation 
effects. Our results show that if one relies on approximate approaches to estimate the 
FSI effects one might underestimate considerably the momentum distributions at high p m 
extracted from experiment in those kinematical regions. 

As stated above, the aim of the present work is mainly to study relative effects of anti- 
symmetrization and FSI, which are often treated approximately, via a complete solution of 
the quantum mechanical few-body problem. We have conducted this study using a semireal- 
istic potential model. Nevertheless it is interesting to compare our results with experiment. 
This comparison is shown in Table 2. Except for the case at the lowest u and q (Kin. N. 
1) where there is a good agreement, our results are almost systematically higher than data. 
The difference ranges from about 30% for the kinematics closer to the quasi elastic ridge to 
about 70 and even 100% for the other ones. This comparison is better illustrated in Fig. 01 
One can see that, while FSI tends to bring theoretical results closer to data for the kine- 
matics at lower momenta (Kin. N. 1-5), it affects in the opposite direction those at higher q 
(Kin. N. 6-9), with the largest effect for Kin. N. 9 which corresponds to the highest q and 
p m -values. This is a delicate region where cross sections are small and potential dependence 
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and relativistic effects neglected here might play a major role. 
B. Non parallel kinematics 

It is interesting to investigate the above effects also in non parallel kinematics. At fixed 
energy and momentum transfer one can access different p m varying 9 p . Therefore in PWIA 
the response reflects a proton momentum distribution. In the following we will show how 
antisymmetrization and FSI can spoil this interpretation. For the (u, q) values of Table |J in 
Fig. 0] results for non parallel kinematics are shown as functions of p m . In the upper panel 
one can clearly see that the mere antisymmetrization effect does not allow the interpretation 
of the response in terms of momentum distribution beyond certain values of p m , depending 
on the kinematics. These values are rather small (around 1 fm _1 ) for the kinematics at 
lower momentum transfer and can reach 2 fm _1 for those at higher q. This is of course 
discouraging for a study of the short range correlations, which contribute mainly to the 
higher tail of the momentum distribution. 

Fig. |3] shows that antisymmetrization effects tend to fill the minimum of the response in 
PWIA. In order to illustrate the FSI effect, in the lower panel we have chosen Kin. N. 3 
with a smaller and Kin. N. 9 with a larger g-value. As in parallel kinematics FSI tends to 
decrease the response in the former case and to enhance it in the latter. It is interesting to 
see that some minima reappear and some are filled when FSI is included. 

For a better understanding of the situation it is instructive to plot the matrix elements 
calculated from Eqs. (jlOjl. (JUjl and ()15|) . As an example we choose Kin. N. 3. In Fig. our 
results for T^ t orn are shown for PWIA and PWIAS. Moreover, in order to see the difference 
between an independent particle model and a correlated one, we also show the corresponding 
results obtained in an harmonic oscillator (h.o.) model. (The h.o. parameters have been 
fixed to the radii of 4 He and 3 H). Since the MTI-III potential has a rather strongly repulsive 
core the comparison exhibits the effect of ground state short range correlations. One readily 
sees that at low p m the MTI-III potential gives a 15% quenching. The tail region is amplified 
in the inset. The results of the two models have similar behaviors with increasing p m , both 
in PWIA and PWIAS (see also inset of Fig. EH)- However, while the h.o. PWIA matrix 
element remains always positive the corresponding one for MTI-III crosses the zero axis, 
giving origin to the minimum visible in Fig. 0] The minimum is then washed out by the 
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antisymmetrization effect. 

In Fig.^Jo the additional role of FSI is shown. In this case the total matrix element T^° rn + 
Tpf 1 is complex. Real and imaginary parts are shown and compared to the Born result with 
the MTI-III potential. In the inset the complicate interplay of the different contributions 
is illustrated. It is evident that FSI leads to a result close to zero for a rather wide p m 
range, causing appearances and disappearances of minima in the cross section. A more 
realistic interaction may change the present picture in that kinematical region considerably. 
Nonetheless this model study points out that it might be difficult to search for ground state 
correlation effects at high p m values within a PWIA picture. 



IV. CONCLUSIONS 

We have presented the results of an ab initio calculation of the 4 He(e, e'p) 3 H longitudinal 
response obtained by means of the integral transform method with a Lorentz kernel. As 
NN interaction the MTI-III potential model is used. The aim has been to investigate the 
limits of the PWIA approximation (factorization in terms of momentum distribution) due to 
the effects of antisymmentrization and FSI. We have analyzed the situation for the parallel 
kinematics investigated in the experiments of Ref. and for two non parallel kinematics. 
Our model study has shown that the factorized approach (PWIA) might be a reasonable 
approximation for small missing momenta (below 0.5 fm _1 ) and higher momentum transfers 
(above 2 fm _1 ). Unfortunately the situation for higher missing momenta becomes much more 
involved. Both antisymmetrization effects and FSI play an important role. In particular for 
non parallel kinematics their entanglement can give rise to drastic deviations from the PWIA 
result. Furthermore, one may expect considerable sensitivity to nuclear dynamics here. On 
the one hand this result can be considered discouraging in relation to the possibility to 
"measure" directly short range ground state correlations. On the other hand it is possible 
that, due to the sensitivity of the response to all effects, those kinematical regions are 
ideal to study potential model dependence, including perhaps that due to thre-body forces. 
However, FSI has to be treated in a proper way and realistic interactions have to be used 
before definite conclusions can be drawn. The integral transform approach with a Lorentz 
kernel is a promising approach to pursue such studies. 
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PWIA 



Kin. 


Q 




Pm 


Fl/(G p e ) 2 


Apwias 


Afull 


N. 


(MeV/c) 


(MeV) 


(MeV/c) 


[(GeV/c) ^sr" 1 ] 


(%) 


(%) 


1 


299 


57.78 


+ 30 


185.2 


+ 9.3 


- 39.6 


2 


380 


83.13 


+ 30 


185.2 


+ 1.2 


- 20.1 


3 


421 


98.19 


+ 30 


185.2 


+ 0.0 


- 12.8 


4 


299 


98.70 


- 90 


100.0 


+ 4.5 


- 43.4 


5 


380 


65.06 


+ 90 


100.0 


+ 3.9 


- 16.6 


6 


544 


126.6 


+ 90 


100.0 


- 1.1 


+ 11.4 


7 


572 


137.82 


+ 90 


100.0 


- 1.9 


+ 11.5 


8 


650 


175.67 


+ 90 


100.0 


- 1.7 


+ 10.8 


9 


680 


146.48 


+ 190 


14.63 


- 4.4 


+ 52.1 



TABLE I: The 4 He(e, e'p) s H longitudinal response function of Eq. © in PWIA approxima- 
tion for the parallel kinematics of Fig. 2. The relative PWIAS and FULL effects, Ax = 
(X - PWIA) /PWIA, are also listed. 
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Kin. F L [(GeV/c)"^!- 1 ] 



No. 




Expt. 




FULL 


1 


59.0 


±2.0 


±2.2 


66.4 


2 


49.6 


±2.1 


±2.1 


66.9 


3 


46.2 


±2.5 


±2.2 


62.7 


4 


27.8 


±1.0 


±1.2 


34.9 


5 


28.4 


±1.2 


±1.3 


37.2 


6 


14.8 


±1.4 


±1.2 


25.9 


7 


16.0 


±1.5 


±1.3 


23.0 


8 


9.96 ±1.29 ±1.15 


16.3 


9 


1.35 ±0.22 ±0.22 


2.73 



TABLE II: The 4 He(e, e'p) 3 R longitudinal response function. Theoretical results (FULL) compared 
to the experimental values of Ref. [17]: statistical and systematic uncertainties are indicated (±stat. 
±syst.) 
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FIG. 1: T = 1, L = [panel aj] and T = 0, L = 2 [panel &J] components of the inclusive response 
(solid line) compared to those of R^f + R n ^ h [see Eq. (JSJ)] (dashed line) for q = 300 MeV/c. 
The arrows indicate the proton-neutron-deuteron (p,n,d,) and deuteron-deuteron (d,d) break-up 
thresholds. 
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FIG. 2: Position of the various kinematics of Tabled with respect to the to = q 2 /(2m p ) ridge. The 
shaded area represents the region below the break-up threshold, where the cross section is zero. 
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FIG. 3: Percentage deviation from the experimental values of Ref. [l7|: PWIA (open circles), 
FULL results (full circles). 
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FIG. 4: The 4 He(e, e'pfR longitudinal response function of Eq. © as function of p m : a) PWIAS 
for the different (uj, q) values of Table [I] labelled with the corresponding numbers; b) PWIAS (dotted 
and dot-dot-dashed line) and FULL (dashed and dot-dashed line) results for the (uj, q) values of 
Kin. N. 3 and 9 of Table [I] The solid line represents the PWIA. 
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FIG. 5: Matrix element T(E) of Eq.@ in non parallel kinematics, for q and u of Kin. N. 3, as 
function of p m . a): h.o. model: PWIA (dashed line), PWIAS (dotted line); MTI-III potential 
model: PWIA (full line), MTI-III, PWIAS (dot-dashed line), b) MTI-III potential model: real 
(dashed line) and imaginary (dot-dashed line) parts of the FULL matrix element. PWIA and 
PWIAS results as in a). 
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